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From an array of quantum wires to three-dimensional fractional topological insulators 
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The coupled-wires approach has been shown to be useful in describing two-dimensional strongly 
interacting topological phases. In this manuscript we extend this approach to three-dimensions, and 
construct a model for a fractional strong topological insulator. This topologically ordered phase has 
an exotic gapless state on the surface, called a fractional Dirac liquid, which cannot be described 
by the Dirac theory of free fermions. Like in non-interacting strong topological insulators, the 
surface is protected by the presence of time-reversal symmetry and charge conservation. We show 
that upon breaking these symmetries, the gapped fractional Dirac liquid presents unique features. 

In particular, the gapped phase that results from breaking time-reversal symmetry has a halved 
fractional Hall conductance of the form axy = if filling is u = 1/m. On the other hand, 

if the surface is gapped by proximity coupling to an s-wave superconductor, we end up with an 
exotic topological superconductor. To reveal the topological nature of this superconducting phase, 
we partition the surface into two regions: one with broken time-reversal symmetry and another 
coupled to a superconductor. We find a fractional Major ana mode, which cannot be described by 
a free Majorana theory, on the boundary between the two regions. The density of states associated 


with tunneling into this one-dimensional channel 
the corresponding Laughlin state. 

PACS numbers: 73.43.-f,05.30.Pr,03.65.Vf, 71.27.+a 


I. INTRODUCTION 

The theoretical predictions and consequent experi¬ 
mental discoveries of topological insulators in two and 
three dimensions [1-12] have enriched our understand¬ 
ing of topological phases of matter. Importantly, it was 
demonstrated that topological states may be much more 
widespread than previously believed, and that in partic¬ 
ular they can exist beyond the realms of two-dimensions 
(2D). Since then, many systems realizing various topolog¬ 
ical phases have been proposed, and a periodic table for 
topological phases of non-interacting fermions has been 
established [13, 14]. While this classification is limited to 
topological states protected by time-reversal symmetry, 
particle-hole symmetry, and chiral symmetry, the exist¬ 
ing tools can be applied to other protecting symmetries. 
A particularly notable extension of the above classifica¬ 
tion are the topological crystalline insulators [15], pro¬ 
tected by the crystal point group symmetries. 

The effects of strong interactions on topological ma¬ 
terials, on the other hand, are much more subtle. Par¬ 
ticularly interesting situations occur when interactions 
stabilize topologically ordered phases that cannot be re¬ 
alized in free-fermion systems. The excitations in these 
systems have fractional statistics, and usually carry frac¬ 
tional quantum numbers. Thus, topologically ordered 
systems are sometimes called fractional phases. The most 
famous example of such a phase is the fractional quantum 
Hall effect (FQHE). 

While topological phases have conclusively been found 
to exist in three-dimensions (3D), fractional phases are 
still usually associated with 2D systems. This is a con- 


proportional to uj^ ^, in analogy to the edge of 



FIG. 1. A schematic depiction of the array of wires studied 
in this manuscript. The coupled-wires approach allows us to 
treat interacting terms using the Abelian bosonization frame¬ 
work, making it very useful in constructing strongly interact¬ 
ing topological phases. Throughout this work, we focus on 
the construction of a fractional strong topological insulator. 


sequence of the well-known theorem stating that anyonic 
statistics between two point-particles cannot occur in 3D 
[16, 17], 

A possible way to go around this theorem and realize 
topologically ordered phases in 3D is to consider loop 
excitations, which can have non-trivial braiding statistics 
with point-particles, as well as other loops. 

Indeed, a few recent works [18-26] have used various 
non-perturbative approaches to discuss fractional topo¬ 
logical insulators in 3D. These are the fractional counter¬ 
parts of the well-studied non-interacting strong topolog¬ 
ical insulators. 

In this work we construct a model for such systems 
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FIG. 2. The fractional strong topological insulator phase we construct has an exotic fractional Dirac liquid on the surface. 
As in the non-interacting case, this gapless surface is protected by time reversal symmetry and charge conservation. Upon 
breaking any of these symmetries, we end up with a gapped surface displaying unique topological properties. For example, by 
breaking time-reversal symmetry, we get a halved fractional quantum Hall effect, characterized by a surface Hall conductance 
of the form axy = where m is an odd integer (in a state associated with the 2D Laughlin state at filling u = 1/m). 

On the magnetic domain wall shown in (a), one therefore finds a chiral Luttinger liquid, similar to the edge mode of the 
corresponding 2D Laughlin-state. If the system is gapped by proximity to an s-wave superconductor (i.e., by breaking charge 
conservation), we get an exotic time-reversal invariant topological superconductor. It is interesting to study the boundary 
between a superconducting region and a magnetic region, as depicted in Fig (b). A fractional Majorana mode, which cannot 
be described by a free Majorana theory, is found on the interface between the two regions. The tunneling density of states 
associated with this ID channel is proportional to , in analogy to the edge of the corresponding Laughlin-state. 


from an array of weakly coupled wires, as illustrated in 
Fig. (1). The ability to introduce interactions, and treat 
them using the bosonization framework has made this 
approach successful in producing fractional phases in 2D 
[27-38]. In addition to reproducing well-known FQHE 
states in the extremely anisotropic limit, it has been 
shown to be useful in producing new fractional states, 
which would otherwise be very hard to realize in terms 
of a microscopic model. 

Here we show that the wires approach can be extended 
to 3D, and use it to construct a fractional strong topo¬ 
logical insulator. 

Since, by assumption, different wires are weakly cou¬ 
pled, we linearize their energy spectra around the Fermi- 
points, and treat inter-wire terms as perturbations within 
this linearized framework. Once the ID spectra are lin¬ 
earized, we can fully describe the single and many body 
excitations using chiral bosonic fields. The only remnants 
of the original model are encoded in the values of the 
Fermi-momenta associated with the various wires (these 
are determined by the overall density, the spin-orbit cou¬ 
pling, the magnetic field, etc.). The Fermi-momenta 
determine, due to momentum conservation, the allowed 
inter-wire terms. To keep track of these terms, it proves 
useful to depict the Fermi-momenta diagrammatically. In 
2D models, for example, we plot the Fermi-momenta as 
a function of an index enumerating the wires (see Figs. 
(9a)-(9b)). Importantly, the diagrams present the prob¬ 
lem as a lattice model in a lower dimension. The task of 
calculating some topological properties of the full model 
is then reduced to the calculation of similar properties in 


the lower-dimensional non-interacting system. For exam¬ 
ple, the analysis of 2D topological insulators is reduced, 
in some aspects, to the analysis of a ID topological sys¬ 
tem. Similarly, the analysis of a 3D strong topological 
insulator is reduced to the analysis of a 2D topological 
insulator. At integer fillings, this enables us to construct 
topological phases in two and three dimensions. 

Upon varying the filling to a given fractional value, 
the pattern formed by the Fermi-momenta is modified, 
and we are forced to consider a different set of inter-wire 
terms, which in general need to involve multi-electron 
processes. At a specific set of filling factors, we find it 
useful to define new fields, associated with a new set of 
effective Fermi-momenta. The transformation to the new 
fields can be chosen such that the corresponding diagram 
is mapped onto the one describing a system at integer 
filling. Then, by repeating the analysis that led to the 
creation of the non-interacting topological phase in terms 
of the transformed fields, we get its fractional analog. 

The resulting 3D topologically ordered phase stud¬ 
ied here will be shown to have a novel gapless surface 
mode, which cannot be described by Dirac’s theory of 
free fermions. Throughout this work, we refer to the sur¬ 
face as a fractional Dirac liquid. In analogy to a strong 
topological insulator, the surface is protected by time re¬ 
versal symmetry and charge conservation. 

It is insightful to study what happens to the fractional 
Dirac theory once it is gapped by breaking any of these 
symmetries. We will see that if a time reversal breaking 
perturbation is introduced, the surface acquires a Hall 
conductance of the form a^y = ^ ^ in a state at filling 
u = 1/m. This Hall conductance is half that of the as- 
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sociated 2D Laughlin FQH state, and we therefore refer 
to this effect as a halved fractional quantum Hall effect. 
This should be compared with the half-integer quantum 
Hall effect that occurs on the surface of a strong topolog¬ 
ical insulator, which corresponds to m = 1 in our frame¬ 
work. As a result of the above, a magnetic domain wall of 
the type depicted in Fig. (2a) contains a gapless channel 
described by the chiral Luttinger-liquid theory, similar to 
the edge states of a Laughlin-state. 

In addition, we will break charge conservation by prox¬ 
imity coupling the surface to an 5-wave superconduc¬ 
tor. The resulting superconducting state is found to 
be topologically non-trivial. On the surface of a non¬ 
interacting strong topological insulator, for example, one 
finds a phase that resembles a + ipy superconductor, 
but has time reversal symmetry [39]. One way to reveal 
the topological nature of the surface is to separate the 
surface into two domains: one with broken time-reversal 
symmetry, and another with broken charge conservation. 
It was shown in Ref. [39] that a chiral Major ana mode is 
localized near the boundary separating the two regions. 
Furthermore, it was shown in Ref. [40] that this remains 
correct in the presence of interactions. 

Repeating the same thought experiment in our frac¬ 
tional phase, we find a chiral self-Hermitian mode on the 
boundary, which cannot be described by a free Major ana 
theory, as illustrated in Fig. (2b). In particular, the 
tunneling density of states associated with this ID chan¬ 
nel is proportional to ^ as opposed to the constant 
density of states characterizing free Major ana fermions. 
Throughout this paper, we refer to this mode as a frac¬ 
tional Majorana mode. 

The structure of the paper is as follows: In Sec. HA 
we discuss the physics of fractional weak topological in¬ 
sulators. These are states constructed by stacking 2D 
fractional topological insulators [41]. While these sys¬ 
tems do not host anyonic excitations which are decon- 
fined in the stacking direction, they possess interesting 
surface states, intimately related to those of the frac¬ 
tional strong topological insulator. Then, in Sec. HI, we 
use an oversimplified yet intuitive model (with a modi¬ 
fied time-reversal symmetry), derived from the surface of 
the fractional weak topological insulator, to describe the 
surface of the fractional strong topological insulator. The 
purpose of this discussion is to pedagogically derive the 
properties of the surface once it is gapped by breaking its 
protecting symmetries. We stress that these properties 
will be derived rigorously in Section IV without using the 
above oversimplified model. 

In the sections that follow, we use a wire construction 
to create a fractional strong topological insulator. The 
analysis is done in a few stages. First, in Sec. IV A we 
review the construction of a 2D fractional topological in¬ 
sulator presented in Ref. [31]. This model will be the 
starting point in our constructions of 3D states. In this 
section we also introduce the general approach and no¬ 
tations used throughout the rest of this work. Then, in 


Sec. IVB-IV C we will show that by stacking copies of 
this 2D phase with a set of non-trivial inter-layer cou¬ 
pling terms, a fractional strong topological insulator can 
be stabilized. Finally, in Sec. IV D, we will use the wires 
model to derive the properties of the resulting phase. 


II. FRACTIONAL WEAK TOPOLOGICAL 

INSULATORS 

A. Weak topological insulators 

Before discussing fractional weak topological insula¬ 
tors, we briefly review the physics of non-interacting weak 
topological insulators. 

To construct a simple model for a weak topological 
insulator, we imagine stacking many 2D topological in¬ 
sulators (see Fig. (3a) for illustration). 

For simplicity, we assume that the 2D topological insu¬ 
lators in the various layers conserve Sz • We can therefore 
describe their edge states as counter propagating spin-up 
and spin-down modes. In the limit where the layers are 
decoupled, the surface in the xz direction, for example, is 
composed of these helical modes. The modes are repre¬ 
sented diagrammatically in Fig. (3b), where the vertical 
direction represents the layer index, and the horizontal 
direction represents the spin. The symbol (0) corre¬ 
sponds to a right (left) mover. 

We now introduce coupling between adjacent layers. In 
doing so, we consider only the terms represented by the 
arrows in Fig. (3b), as these are the only terms capable of 
gapping the surface. Time reversal symmetry relates the 
amplitudes of the terms represented by the full arrows 
and those represented by the dashed arrows. 

To be concrete, if we define the electron annihilation 
operators ' 0 n,t(^) 5 ' 0 n 4 (^) (where n is the layer index), 
we can write the low energy surface Hamiltonian in the 
form 

H = H, + H,, (1) 

where 

Hx = -iv'Y^ / dx (2) 

represents the Hamiltonian of the decoupled helical edge 
modes, and 

j dx {^l,^{x)i)n+i,i{x)e^°‘ (3) 

n 

+ h.c)j 

represents their nearest-layer coupling. In the above, a 
is a fixed phase, which we set to be 0 for simplicity. 
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FIG. 3. (a) The simplified model we use to describe a weak topological insulator. Each plane forms an conserving 2D 
topological insulator. The edge of each layer therefore contains counter-propagating spin-up and spin-down modes, represented 
by red and blue arrows, (b) A diagrammatic depiction of the surface of the simplified model. The vertical direction represents 
the layer index, and the horizontal direction represents spin. The symbol (g) (Q) corresponds to a right (left) mover. We 
introduce nearest-layer coupling terms, represented by arrows connecting different chiral modes. Time reversal symmetry relates 
the amplitudes of the terms represented by the full arrows to the amplitudes of the terms represented by the dashed arrows. 
The above model can be decomposed into two decoupled gapless Hamiltonians, represented by the green and purple arrows. 


If the system has periodic boundary conditions in the 
X and ^ directions, kx and are good quantum numbers. 
We can therefore write the Hamiltonian in /c-space, where 
it takes the form H = J dkx^p^ {k)'H{k)^{k), with 
ip{k) = (V’4,(fc) and 

H = —2t(Ty sin {kzo) + vkx<Jz (4) 

(here the a^’s are Pauli matrices acting on the spin de¬ 
grees of freedom). We see that we have two Dirac cones 
on the surface: one around {kx = O^kz = 0), and another 
around {kx = O^kz = tt). These are not protected by 
time-reversal symmetry, as time-reversal invariant terms 
can couple the two Dirac modes and gap them out. How¬ 
ever, such terms necessarily involve a large momentum 
transfer. Therefore, if lattice translation invariance is 
also imposed, the two Dirac modes remain protected. 

B. The surface of fractional weak topological 
insulators 

In order to generalize the above simplified model to a 
fractional weak topological insulator, we imagine stack¬ 


ing layers of 2D fractional topological insulators. Again, 
we focus on the Sz conserving case, where each 2D frac¬ 
tional topological insulator can be thought of as two de¬ 
coupled FQHE layers. These two FQH states have equal 
densities, but opposite fillings ±iy and spins. 


For now, we keep the discussion general and do not 
provide a specific model for the 2D fractional topologi¬ 
cal insulators in the various layers. For such a concrete 
model, the reader is referred to Sec. IV A, where we 

review the wire construction of a fractional topological 
insulator introduced in Ref. [31]. 


We focus on the time reversal invariant analog of a 
Laughlin state with filling u = 1/m, where m is an odd 
integer. Before introducing inter-layer coupling terms, 
the natural way to describe the edge modes of the various 
layers is in terms of two boson fields Xn Xn- These 
modes satisfy the commutation relations 


Xn{x),Xn'{x') 


= —i'n:aSn,n'Sa,a'sigrL{x - x') + T-j,rsign(n - n') + 2—TT5n,n'(^Z'‘’ 


( 5 ) 


where in the above a = 1(—1) for spin t (i)- terms of these, the chiral fermion operators take the form 
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'0n,cr oc and the Laughlin quasiparticle operators 

are given by ^ If we introduce inter-layer 

coupling terms, the effective low energy Hamiltonian de¬ 
scribing the surface is given by 

H = H, + H,, (6) 

where 

( 7 ) 

n,cr 

represents the decoupled chiral Luttinger liquids, and 
Hz = f dx 

n ^ 

- + ^-C-) 

represents their coupling. In what follows we assume that 
the amplitude i is large enough such that these operators 
flow to the strong coupling limit (if they are considered 
separately). In the integer case, where m = H reduces 


to the Hamiltonian H defined in Eq. (1), as can be seen 
directly by rewriting H using Abelian bosonization. In 
the fractional case, where m > H cannot be mapped to 
a non-interacting fermionic Hamiltonian. However, the 
model can still be represented by the diagram shown in 
Fig. (3b). Notice that now the symbols 0 and 0 rep¬ 
resent right and left moving (or y) fields, respectively. 

Each of the two terms in Eq. (8), if considered sepa¬ 
rately, can gap out the spectrum in the thermodynamic 
limit. However, the two gapped phases that result from 
these terms are topologically distinct. Noting that the 
two non-commuting operators have the same amplitude 
and scaling dimension, it is clear that the system must 
be in a critical point between the two gapped phases. We 
emphasize that the criticality of the surface is imposed 
by time reversal symmetry, and is not a result of fine- 
tuning the Hamiltonian. It is explicitly assumed that 
additional interacting terms do not destabilize the criti¬ 
cal point. This crucial assumption must be checked for a 
given microscopic model. 

Like in the integer case, this gapless model is rep¬ 
resented by two decoupled gapless theories. This can 
be identified by writing the Hamiltonian in the form 
H = Hi + ^2, with 


^ E / (^:rXL) + {dxXL+l) + tY j (^L,t^2n+14 - ^L+2,t^2n+14 + h.c)j , (9) 

^2 = ^ E / (^^xL+l) + {dxxL) + tY j (^L+24^2n+l.t - ^L 4 ^ 2 n+l.t + h.c)j . (10) 

n,cr n 


The green (purple) arrows in Fig. (3b) represent terms 
belonging to Hi (^2)- In the integer case each of these 
Hamiltonians is described by a low energy Dirac theory, 
as can be seen by refermionizing the bosonic theories. Im¬ 
portantly, the Hamiltonians Hi and H 2 remain gapless in 
the fractional case as well, as indicated by the argument 
we have used to show that H is gapless. However, in the 
fractional case these Hamiltonians cannot be described 
by a low energy Dirac theory. In the future we refer to 
their low energy theories as fractional Dirac liquids. 


III. A SIMPLIFIED MODEL FOR THE 
SURFACE OF FRACTIONAL STRONG 
TOPOLOGICAL INSULATORS 

We have seen above that the simple model for a frac¬ 
tional weak topological insulator can be decomposed into 
two decoupled gapless theories described by the Hamil¬ 


tonians Hi and H 2 . In the integer case, weak topologi¬ 
cal insulators possess an even number of Dirac cones on 
the surface, while strong topological insulators necessar¬ 
ily have an odd number of Dirac cones [10]. By analogy, 
we naturally expect the low energy theory describing the 
surface of a fractional strong topological insulator to be 
related to the surface of a fractional weak topological 
insulator. We therefore ask whether a single decoupled 
surface Hamiltonian, say iLi, describing left movers on 
layers with even indices and right movers on layers with 
odd indices (see Fig. (4a)), can faithfully describe the 
surface of a fractional strong topological insulator. 

Strictly speaking, this is clearly impossible since Hi is 
not independently time reversal invariant. Nevertheless, 
we follow Ref. [34] in noting that Hi is invariant under a 
modified time reversal operation, defined as the product 
of the original time-reversal operator and a translation 
by a unit cell. 

Such a symmetry characterizes antiferromagnetic 
topological insulators, which can be created, for exam- 
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FIG. 4. (a) The simplified model we use to describe an antiferromagnetic topological insulator and its fractional analog. Each 
layer can be thought of as a single-component QHE system at Filing (where z/ = 1 in the integer case, and in the fractional 
case we consider v — Ijm^ where m is an odd integer). Due to the existence of a modified time-reversal symmetry, given by 
the product of the original time-reversal operator and a translation by a unit cell, the xz surface remains gapless even in the 
presence of inter-layer coupling terms. The surface can be described in terms of coupled ID channels, which enables a simple 
derivation of universal surface properties. We argue that these properties remain true for a strong topological insulator and its 
fractional analog, where the surface is protected by the local time reversal symmetry, (b) A diagrammatic representation of xz 
surface of the above model in terms of the coupled ID chiral edge modes, (c) The diagram that corresponds to the situation 
where the modified time-reversal symmetry is explicitly broken by introducing a term of the form (11) with t' = t. 


pie, by adding an antiferromagnetic order parameter to 
a strong topological insulator (without closing the gap). 
Remarkably, it was shown in Refs. [42, 43] that the in¬ 
troduction of such a time-reversal breaking perturbation 
does not destroy all the surface Dirac cones. Instead, 
the remaining Dirac cones are protected by the modified 
time-reversal operation described above. 

This leads us to assume that Hi faithfully describes 
the surface of the fractional analog of an antiferromag¬ 
netic topological insulator. The corresponding gapless 
surface is expected to be in the same universality class 
as the surface of a fractional strong topological insulator 
with a local time-reversal symmetry. We can therefore 
use Hi to derive some of the universal properties ex¬ 
pected to characterize the surface of a fractional strong 
topological insulator. Notice that the bulk excitations of 
the fractional antiferromagnetic topological insulator are 
generally different from the fractional excitations charac¬ 
terizing the 3D fractional strong topological insulator. 

As noted in Ref. [34] , studying a system with the mod¬ 
ified time-reversal symmetry, for which one can write an 
explicit model of the surface in terms of weakly coupled 
ID channels, greatly simplifies the analysis. This is im¬ 
possible in systems that have a local (unmodified) time- 
reversal symmetry. To study these using a set of weakly 
coupled ID systems, we will be forced to explicitly con¬ 
struct the 3D bulk as well. This will be done in Sec. 
IV. 

Again, we depict the surface model using diagrams rep¬ 
resenting the inter-layer terms coupling the various chi¬ 
ral modes. In particular, the diagram that represents the 


gapless surface Hamiltonian Hi is shown in Fig. (4b). 

Until now we have preserved the modified time-reversal 
symmetry and charge conservation of Hi^ and thus the 
surface was found to be gapless. In what follows we focus 
on the properties of the surface once these protecting 
symmetries are explicitly broken. We will see below that 
the resulting gapped fractional Dirac mode has unique 
properties. 

To break the modified time reversal symmetry, we add 
a perturbation of the form 


n 

+ '0L+2,t'^2n+l,t + (11) 

The physics of the resulting phase becomes transparent 
at the points i' = where the Hamiltonian is decom¬ 
posed into decoupled sine-Gordon models. For example, 
if t' = t the non-quadratic part of the Hamiltonian takes 
the form 


V / + h.c)j = 


At' / dx cos m ^ 


^ U2n -X2n+1 


( 12 ) 


Thus we end up with a set of mutually commuting cosine 
terms, which can be pinned in the ground state and fully 
gap the surface. This configuration is depicted diagram- 
matically in Fig. (4c). 
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FIG. 5. A diagrammatic representation, in terms of the (or x) modes defined at the beginning of Sec. IIB, of the situation 
where the coefficient i' of the time reversal breaking term changes abruptly from —i to i around n = 1. A decoupled chiral 
-0 mode is seen to be localized near the boundary at n = 1, leading to the conclusion that the gapped surface has a halved 
fractional quantum Hall effect. 


To identify the surface Hall conductance axy charac¬ 
terizing this gapped state, we change i' abruptly from 
—t to t around the point n = 1, as demonstrated in Fig. 
(5). It can be seen that a chiral x-mode, localized on the 
boundary, remains decoupled. This mode is identical to 
the mode residing on the edge of a z/ = 1/m Laughlin 
FQHE state. By invoking the bulk-edge correspondence 
and noticing that this mode is a result of contributions 
from the two sides of the boundary, we find that the 
gapped surface has axy = We note that since 

the perturbation defined in Eq. (11) is small compared 
to the bulk energy gap, the resulting gapless chiral mode 
and the associated value of axy must be interpreted as 
properties of the surface. 

Thus, the gapped surface of the fractional strong topo¬ 
logical insulator exhibits a halved fractional quantum 
Hall effect. We emphasize that the above result is not 
limited to the special choice of parameters, which were 
tuned to the so called “sweet point” [44, 45], and is in 
fact true as long as the gap remains open. Physically, 
deviations from f in the two regions introduce a 

non-zero localization length for the boundary mode. 

Alternatively, the fractional Dirac liquid can be gapped 
out by proximity-coupling to an s-wave superconductor, 
i.e., by breaking charge conservation. A superconducting 
term which does not violate the modified time-reversal 
symmetry is given by 


Ha 



'02n,t'^2n+14 


+ '02n,t'02n-l,4, + h.C. 


(13) 


Again, the physics becomes simple at a sweet point 
given by A = Taking A = t, for example, the non¬ 
quadratic part of the Hamiltonian can be written as 

Ha = (- 1 )” f ( 14 ) 


with ipn,a = 7n + ( 7 ^ ^ (^otice that we 

have omitted the spin index, which is fully determined by 
n). The modes 7 ^ are referred to as fractional Majorana 
modes. Indeed, in the integer case, these become chiral 


Majorana modes, described by a free Majorana theory. 
The structure of this Hamiltonian is depicted diagram- 
matically in Eig. ( 6 b), where the dotted symbols repre¬ 
sent the fractional Majorana modes 7 ^ and 7 ^. Erom the 
form of Eq. (14), it is clear that the system is gapped. 
Notice that Eq. ( 11 ), describing time-reversal symmetry 
breaking, can be expressed in terms of the fractional Ma¬ 
jorana modes as well. In particular, the non-quadratic 
part of the Hamiltonian Hi + Hf with i' = i takes the 
simple form 2ii / dx ( 72 n+i 72 n “ 72n+i72n) > which is 
depicted in Eig. ( 6 a). 

Once the surface is gapped by proximity to a supercon¬ 
ductor, it forms an exotic time-reversal invariant topolog¬ 
ical superconductor. We leave the further investigation 
of such a novel superconducting state to future works. 

However, it is illuminating to examine the boundary 
between a magnetic region with broken time-reversal 
symmetry, and a superconducting region. The physics 
is simplest if we gap the region with n > 1 using a time 
reversal breaking term of the form ( 11 ) with i' = t, and 
the region with n < 1 using a proximity-coupling term 
of the form (13) with A = i. This situation is depicted 
in Eig. ( 6 c). We see that a decoupled chiral Majorana 

mode of the form 7 i = | is localized 

on the boundary. The propagator describing this frac¬ 
tional Majorana field takes the form (7K^5 ^)7i (O 5 0)) ^ 
{x-\-vt)~^, making it clear that this gapless channel can¬ 
not be described by a free Majorana theory. We point 
out the similarity to the edge mode of the 2D fractional 
topological superconductor discussed in Ref. [46]. 

One can generalize the construction to study hierar¬ 
chical states, with u different from 1/m. Within the 
simplified model, each layer in Eig. (4a) now contains 
a EQHE state with a general iF-matrix of dimension d, 
and therefore d distinct edge modes Xn (where n is the 
layer index and j = 1, • • • , d enumerates the edge modes 
in a given layer). Switching on terms that preserve the 
modified time reversal symmetry and couple only modes 
with the same j, we get a generalized fractional Dirac 
theory associated with any iF-matrix. It is clear, how¬ 
ever, that the generalized fractional Dirac theories are 
not necessarily protected by symmetries, as some modes 
can be gapped out by coupling fields with different j’s 
(without breaking time-reversal symmetry). 

If we repeat the analysis presented in this section and 




break the modified time-reversal symmetry (still with¬ 
out coupling modes with different j’s), we get a halved 
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fractional Hall conductance of the form (Jxy = ^ • Addi¬ 
tional coupling between modes with different j’s has the 
potential of changing the Hall conductance. Therefore, if 
symmetry breaking perturbations are introduced, a num¬ 
ber of topologically distinct gapped phases may arise on 
the surface. 

The above analysis relied on a modified time reversal 
symmetry to directly model the surface using coupled 
ID modes. In what follows we treat a 3D model with a 
local time reversal symmetry. As we will see, the analysis 
presented in the next sections reproduces the universal 
results derived here. 


IV. FRACTIONAL STRONG TOPOLOGICAL 
INSULATORS 

A. Construction of 2D fractional topological 
insulators 

In Sec. IVB-IVC we will construct a 3D model for 
a fractional strong topological insulator. Our starting 
point will be the 2D construction of a fractional topolog¬ 
ical insulator which was introduced in Ref. [31]. In this 
section we review this construction, and introduce the 
general approach used throughout the rest of this work. 
Unlike the pervious section (Sec. HI), where we used a 
simplified model with a modified time-reversal symmetry 
[see Fig. (4a)], here we have a local (unmodified) time 
reversal symmetry. To avoid confusion, we will use a 
different set of notations in the analysis that follows. 


1. Two-dimensional topological insulator from an array of 
quantum wires 

The 2D model which is the focus of this section is made 
of 4V wires, as depicted in Fig. (7). The j’s unit cell 
is composed of four wires. We note that for convenience 
the unit cell is shifted by one wire with respect to the 
unit cell defined in Ref. [31]. 

The first and last wires in each unit cell contain elec¬ 
trons, i.e., their highest occupied states are close to the 
minimum of the conduction band. In a particle-hole sym¬ 
metric fashion, the two other wires contain holes, i.e., 
their highest occupied states are close to the maximum 
of the band. 

We describe the system in terms of a tight binding 
model, where each wire is composed of sites at positions 
X = a • i (see Fig. (7)). Here a is the distance between 
adjacent sites, and i is an integer enumerating the sites. 
Adjacent lattice points are coupled with a hopping am¬ 
plitude tx {—tx) in the electron (hole) wires, as depicted 
by the arrows in Fig. (7). We define the annihilation 


operator for an electron with spin s in wire num¬ 

ber a {a = 1.. A) of the unit cell labeled by (i, j). The 
Hamiltonian that describes the decoupled wires takes the 
form Hq = Hm + Hx , where 



and 

= -4 E f E (16) 

i,j \a=l,4 

- E + h.c. 

a=2,3 


The term Hx describes intra-wire hopping, and the 
term Hm produces an opposite shift in energy for the elec¬ 
trons and the holes. We note that to have a particle-hole 
symmetry, the chemical potential is set to be zero. Going 
to /c-space, we write H = with 

i) 2 t V’3t ^4t V’ 24 , i’ai In addi- 

tion, we define the matrices (i = 1... 3) as the Pauli- 
matrices operating on the wires 1-2 (and 3-4) space. Sim¬ 
ilarly, the matrices operate on the (l,2)-(3,4) blocks, 
and Si act on the spin degrees of freedom. 


In terms of these, the Bloch Hamiltonian correspond¬ 
ing to Eqs. (15)-(16) takes the form 

Ho = [m- 2tx cos {kxa)] (17) 


Furthermore, we introduce Rashba spin-orbit interac¬ 
tions with an alternating coupling Ago (~1)^- If the elec¬ 
tric field is aligned in the y direction, the resulting term 
is 

Hso = ^so^z'^z sin (^kxCi) . (I^) 

It proves convenient to define new parameters: 
tx , kso , and k ^, satisfying 

tx = tx cos (ksoa ), (19) 

Ago = ‘2tx sin (ksoa ), 
m = 2tx cos [k^a) . 

In terms of these, it is simple to see that the energy 
spectra of the four decoupled wires in each unit cell take 
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FIG. 6 . A diagrammatic representation of the surface model in terms of the fractional Majorana modes 7 ^ (defined after Eq. 
(14)) in three situations: (a) The modified time reversal symmetry is broken by the perturbation ( 11 ) with ? = i. (b) Charge 
conservation is broken by a superconducting term of the form (13) with A = t. (c) The region with n < 1 of the surface is 
gapped by a superconducting term, and the region with n > 1 is gapped by a modified time-reversal symmetry breaking term. 
A chiral fractional Majorana mode is attached to the boundary at n = 1. Here the vertical direction represents the index z, 
and the horizontal direction represents the layer index n. Dotted symbols represent the chiral fractional Majorana modes. 


the form 

Ei^s = 24 [cos {kPpa) - cos - sk^o) a)] , 

E2,s = -2G [cos (kpo) - cos {{ka: - skso) a)] , 

Es,s = -24 [cos {k%a) - cos ((/c^ + sk^o) a)] , 

^4,s = 2G [cos [kpa) — cos {{kx + skso) cl)] • (20) 

If we define the filling factor as 


picted in Fig. (8a) for the u = 1 case, and in Fig. (8b) 
for the = 1/3 case. We now introduce small tunnel¬ 
ing operators that couple adjacent wires, and write the 
inter-wire Hamiltonian in the form Hy + Hy, where 


H, = -t, E vxf ++ h.c. 


( 22 ) 

K = -<; E + I'.c] , (23) 


hJ 



( 21 ) 


the spectra corresponding to the spin-up sector are de- 


with ty^ty tx. The parameter ty describes hopping 
between two electron-wires or two hole-wires, whereas t'y 
couples the electron and hole wires. 

The Bloch Hamiltonian can now be written in the form 


1-L = 2tx [(cos (/c^a) — cos {kxo) cos {ksoo)) cfzTz — SzTz sin {ksoo) sin {kxa)\ 

- - h (Tyay + TxCf^) - h - TyGy) cos (4fcj,a) - h {j^Oy + TyG^) sin {4kya) . (24) 


For simplicity, we treat the case where the lattice spac- 
ings are identical in the two directions: ax = ay = a. 
This requirement can be lifted without affecting any of 
the topological properties. Notice that in these con¬ 
ventions, the first Brillouin zone is defined as ky G 


3^] ^ (“f ’ f] • investigate the in- 

teger case, = 1, where it can be checked from Eq. (24) 
that as long as ty 7^ ty and 7^ 0, the system is com¬ 
pletely gapped when periodic boundary conditions are 
employed. At ty = ty the gap closes, indicating that there 
might be a phase transition between two topologically 
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FIG. 7. The wire construction we use as a realization of a 
2D fractional topological insulator. This model is the starting 
point in our construction of a 3D fractional strong topological 
insulator (see Sec. IVB-IVC). 


distinct phases. To understand the nature of the insulat¬ 
ing phases for different values of ty and ty, it is useful to 
linearize the spectrum around the Fermi-momenta. 

Assuming that ty^ty tx, the inter-wire terms can 
be treated as perturbations within the linearized the¬ 
ory. Once the spectrum has been linearized, the only 
remnants of the original model are the Fermi momenta 
of the various modes. If momentum conservation is im¬ 
posed, the values of the Fermi momenta severely restrict 
the allowed inter-wire coupling terms for a given u. To 
keep track of these terms, it is convenient to present the 
Fermi-momenta as a function of the wire index (i, a), as 
depicted in Fig. (9a) for the spin-up sector. 

Again, the symbol 0 (0) represents a right (left) 
moving mode, and the arrows represent the coupling be¬ 
tween different modes, generated by the terms defined 
in Eqs. (22)-(23). Since the system is fully gapped 
for 0 < ty < ty, any such state is adiabatically con¬ 
nected, and therefore topologically equivalent, to the 
state where ty is negligibly small compared to ty. There¬ 
fore, it is simple to see that the terms + and 

-^h.c. (and in the same way, 

and + h.c. for the spin-down sector) domi¬ 


nate, and gap the modes near kx = 0. These terms, how¬ 
ever, leave two gapless modes on each edge of the system 
- one for each spin. For the edge at i = 1, these edge 
modes are , and The terms h.c. 

and + (and in the same way t 

and + h.c.) are responsible for gapping the 

modes with kx 0. We therefore have a fully gapped 
bulk with counter-propagating edge modes, protected by 
time reversal symmetry. Thus, the phase with ty < ty is 
a topological insulator. 

On the other hand, any state with ty > ty is adiabati¬ 
cally connected to the state with = 0, where now the 
terms + h.c. and -h h.c. (and in the 

same way '^3 

modes near kx = 0. In this case we have a fully gapped 
bulk with no edge modes, indicating that the phase we 
discuss is topologically trivial. 


2. Two-dimensional fractional topological insulator from an 
array of quantum wires 


Having found that the above model (defined in Eqs. 
(15)-(16) and (22)-(23)) can be a topological insulator in 
the u = 1 case, we now turn to study fractional fillings. 
We restrict ourselves here to fillings of the form u = Ijm, 
where m is an odd integer. More general situations were 
discussed in Ref. [31]. 

The diagram that corresponds to the = 1/3 case is 
shown in Eig. (9b). It is evident that simple tunneling 
processes do not conserve momentum and are therefore 
incapable of gapping the system in this case. Since we 
are interested in a gapped phase, we have to introduce 
interactions and consider multi-electron processes. 

Eortunately, interaction terms become manageable if 
we use the standard Abelian bosonization technique. We 
define the chiral boson fields , such that (x 


exp 


i,RlL 


■ti.RlL 




where k^d^J^ is the corre¬ 


sponding Eermi-momentum, and the boson fields satisfy 
the commutation relations 


/n ,p 

Tolls' 


— ipud 5qi^^q^i 5p p! X ) 

+ m (sign(n - n') + (5„,„/sign(a - a')) + (5„,„'<5„,a'7r (ct®’* + ) , 


(25) 


define new chiral fermion operators according to 




i,R/L ^ 
a,s 




i,R/L 


i,L/Ry 


oc e 




with 


(26) 


i,R/L ^ 
^a,s 


m + 1 


^i,R/L rn-1 ix/R 

'ra,s 2 ^a,s 5 

(27) 

,J,R/L m-'^jJXIR 

^a,s 2 ^ 0 !.,s 

(28) 
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FIG. 8. The spectrum of the system depicted in Fig. (7), when the inter-wire terms are switched off for (a) v — 1, and (b) 
z/ = 1/3. We note that only the spin-up sector is presented here. The energies in blue, cyan, red, and green correspond to wires 
number 1, 2, 3, and 4, respectively. 


We note that Eq. (26) makes sense only when the 
Fermionic operators of the same type are located at close 
but separate points in space. The commutation relations 
of the new bosonic ?7-fields are 




n ,p 


= vm7r(5s,s/(5„,a'<5p,p'<5n,n'sign(a; - x') 

+ in (sign(n - n') + (5„,„/sign(Q; - a')) + (5„,„'(5a,a'7r (cr®’® + . 


(29) 


Plotting the diagrams corresponding to the g’s, we find 
that the picture is identical to the one associated with 
u = 1 (Fig. (9a)). The transformation from the original 
Fermionic degrees of freedom to the composite '0-fields, 
which takes a simple linear form in terms of the boson 
fields, can therefore be interpreted as a transformation 
from z^ = l/mtoz^ = l. Hence, in terms of the '0-fields, 
we can repeat the analysis of the = 1 case, and write 
the terms used to obtain a topological insulator. Writing 


these terms using the r] bosonic fields, we have: 

Ht=J2 j cos (t?!^ - +R^L, 


dxty COS 


( i,R i,L\ 


dxt'y cos 


d 

+ E E / 

+ E E / (%’,s - vZ) 

s p=R,L 


R i —)■ 

(iK-fg) 


(30) 


Notice that while these are simple tunneling opera¬ 
tors in terms of the 0-fields, they describe multi-electron 
processes of the type depicted by the arrows in Fig. 
(9b), in terms of the original 0-fields. It was shown in 
Refs. [27, 28, 47] that there is a large range of strong 
density-density interactions for which such terms flow to 
the strong coupling limit. Furthermore, for all density- 
density interactions, the terms t'y and ty flow to the 
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FIG. 9. The diagrams which present the Fermi-momenta of 
the 2D model as a function of the wire index (labeled by (z, a), 
where i and a are indices enumerating the unit cells and the 
position within the unit cell, respectively) for (a) z/ = 1 , and 
(b) v — 1/3. Notice that only the modes with spin-up are 
presented here. 


strong coupling limit if their bare values are large enough, 
as their flow diagram coincides with that of the Sine- 
Gordon model. 

Similar to the integer case, if ty < ty the first two terms 
in Eq. (30) dominate. Assuming they are made relevant, 
it is evident that the bulk becomes gapped and that the 
two modes and remain decoupled. Each of these 
modes is similar to the edge mode of a z/ = 1/m Laughlin 
QHE state. One can now follow Ref. [28] and show that 
the bulk excitations are fractionally charged and have 
fractional statistics. 

The above analysis suggests that the fractional analog 


of a given integer topological phase can be realized if we 
manage to construct the integer phase from an array of 
coupled wires. Once this is done, we can obtain a frac¬ 
tional phase by considering an analogous system where 
the z/z-fields are replaced by the z/^-fields (or equivalently, 
(j) is replaced by 77 ), as demonstrated above. 

This motivates us to use a similar approach in the con¬ 
struction of fractional phases in 3D. Toward this goal, 
in the next section we will construct a non-interacting 
strong topological insulator from an array of wires. Then, 
in Sec. IV C we will study the fractional phase obtained 
by reducing the filling to v = 1 /m, and constructing a 
strong topological insulator in terms of the composite 
fields. 


B. Strong topological insulators from weakly 
coupled wires 

In order to construct a strong topological insulator we 
stack 2D layers, each made of the wire construction pre¬ 
sented in Sec. IV A, and tune the system such that u = 1. 
The resulting 3D system is made of an array of wires, as 

illustrated in Eig. (1). Eor simplicity, we assume that 
the distance between adjacent layers is a as well. The 
goal of this section is to engineer a time-reversal invari¬ 
ant system with a single Dirac cone near the first and 
last layers. 

To do so, we start by tuning each layer to the critical 
point between the topological and the trivial phase, such 
that the 2D bulk contains two Dirac modes. This can 
be achieved by taking ty = ty, in which case the Bloch 
Hamiltonian describing a single layer is given by 


^xy = 2 G [(cos (ksod) - cos cos (ksod)) CFzTz “ siu (ksod) sin (kxd)] 

(V, + r.a.) - I (r.,. - r,„,) cos (4fe,o) - | (t.o, + T,„.) sin {ik,a ). (31) 


To create a topologically non-trivial gapped 3D phase, we 
perturb the above gapless Hamiltonian by an inter-layer 
term of the form 

^ “ 2G cos {kzd)) Tx (32) 

+2t'^ sin {k^a) SyT^] (1 - a ^), 

with m' ,tz^t'^ ty tx. Unless otherwise noted, we 
implicitly assume that all the coupling constants are pos¬ 
itive. 


a strong topological insulator, we now show that if the 
system is cut at the z = 0 plane, a single Dirac mode is 
localized near the resulting surface. It is clear that as long 
as the inter-layer coupling terms are small, the important 
degrees of freedom are those close to the Dirac points in 
each layer. We therefore project the Hamiltonian onto 
the low-energy subspace of the intra-layer Hamiltonian 
(Eq. (31)). 

The two Dirac cones are located at kx = ky = 0. 
Therefore, to identify this low-energy subspace we solve 
the equation 


To see under which circumstances this model forms 


'kixyikx — ky — O^'lpQ — 0 


(33) 
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for the vectors t/^o- The resulting subspace 
can be spanned by the four vectors |dil,5) = 
\(7x = —l,Ta^ = defined such that |il,«s) = 

-\±l,s) ,Tx\±l,s) = zb |zbl, 5 ) , 5';^ |zbl, s) = 5|zbl,s). 
In what follows, the basis vectors are ordered as 
{|l,t) , I-I 4 ) , |14) As expected, the xy 

Hamiltonian, projected onto the low-energy subspace 


and expanded to first order in momenta, takes the form 
of two Dirac cones. To be specific, in the above basis 
the intra-layer Hamiltonian takes the form 


'Hxy = 


0 h]3irac \ 
^Dirac 0 / 


(34) 


where 


^Dirac — 


0 

2tx sin (ksoa) kxa — 2itykyL 


—2tx sin (ksoa) kxa + 2itykyL 


0 


(35) 


Once the inter-layer part of the Hamiltonian is pro¬ 
jected onto the same low energy subspace, it takes the 
form 




^ f hiD{kz) 0 


V 0 


KM 


where hi-D{kz) = d • a, 
(0, 2t'^ sin (kza ), m' — 2tz cos (kza))^ 
of as the Hamiltonian a ID model. 


can 


with 
be 


(36) 


d = 
thought 


If m' < 2tz , the planar vector d winds once around the 
origin as kz winds around the Brillouin zone, indicating 
that the ID model is topologically non-trivial. Indeed, in 
this regime the model can be shown to have zero-energy 
end modes. 


Since the i part of our 3D Hamiltonian consists of the 
two decoupled models hiD and Kd^ produces two zero- 
energy modes on the z = 0 surface. It is clear that the 
full low energy Hamiltonian, projected onto the subspace 
spanned by these two end modes, must be that of a single 
Dirac cone. 


To see this explicitly, we need an analytic form for 
the two end modes. This becomes simple if we restrict 
ourselves to the regime where 2tz is close to m and the 
gap is small. In this limit, focusing again on the low 
energy properties, we can expand Eq. (36) in small kz, 
and get a continuum model. The model hm then takes 
the form 


Kd 


m' — 2tz — tzd‘1 —2t'^adz 

2t'^adz — ija' — ~ 


(37) 


in the position basis. 


Assuming that the system ends at 2 ; = 0, we look for 
localized zero-energy eigenstates of this Hamiltonian and 
its complex conjugate, which vanish at z = 0. Plugging 
in an exponentially decaying function as an ansatz, we 


find two such solutions 


V2 


/W,|2) = ^ 



f(z), (38) 


with f{z) = \+z _ g and \± = 

—2tz)tz 

atz 


Projecting the xy part of the Hamiltonian (Eq. (34)), 
onto the subspace spanned by the two end modes (Eq. 
(38)), we find a single anisotropic Dirac cone on the sur¬ 
face. The corresponding Hamiltonian takes the form 

'kCxy-suria.ce — iJ^y^y^y ~b kxtx ^1^ {k so^) O'x ] • (39) 


Two remarks are in order: first, one can derive the 
topological nature of the model from the bulk wavefunc- 
tions. In fact, recognizing that the system has an inver¬ 
sion symmetry, generated by the operator (TxTx, we can 
implement the procedure introduced in Ref. [II] (relat¬ 
ing the Z 2 invariant to the parity of the occupied states 
at the time-reversal invariant points) and easily evaluate 
the Z 2 invariant. We have verified that the results of this 
analysis are consistent with the above derivation, where 
the surface was studied directly. 


Second, we note that since the system is guaranteed to 
preserve its topological nature as long the gap remains 
finite, the various parameters are not restricted to the 
values given above. In particular, the strict requirement 
ty = ty can be relaxed. 

In the next section we study the fractional analog of 
the model introduced in this section. 
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C. Fractional strong topological insulators from 
weakly coupled wires 

In Sec. IV A, when constructing the Laughlin-like v = 
1/m state in 2D, we saw that one can define the ry-fields, 
in terms of which the problem is mapped to the simple 
V = \ case. Reversing the logic, we see that by taking a 
topological state with z/ = 1 , and replacing the 0 fields 
by the 77 fields, we expect to get a fractional state. 

In this Section, we follow this approach in generalizing 
our strong topological insulator to its fractional analog. 
To do so, we start by writing the Hamiltonian of the = 1 
case, discussed in the previous section, in terms of the 
bosonic 0-fields. Then we write the same Hamiltonian 


with the 77 -fields, and tune the system to filling z^ = 1 /m, 
where such terms conserve momentum. 

We note that all the fields which are not around = 
0 are trivially gapped by the intra-layer terms, as seen 
most clearly in Fig. (9a). Consequently, the topological 
properties involve only the modes around kx = 0. In 
what follows, we therefore write only the operators acting 
on these fields. In addition, we omit the indices R/T, 
which are fully determined by the spin indices for the 
kx = 0 modes. 

We define the vector 

(40) 

where i represents the index enumerating the unit cells 
in each layer, and n counts the layers. In these notations, 
the low energy Hamiltonian takes the form 




(41) 


H, 


- tyJ2l 


dx 


+ \ + TcoCTx)] + iTy) + IGy) + h.C. 


(42) 


Hz = J dx (^[„+i {-tzTx + it'zSyTz) (1 - (Ta;) ^i,n + h.C^ + (1 - a^) 


(43) 


For simplicity, we do not consider the effects of density- 
density interactions between the various modes. 

We emphasize that analyzing the problem directly in 
terms of the bosons is essential in the fractional case. Un¬ 
fortunately, the bosonic form makes it difficult to see that 
the above set of non-commuting terms results in a gapped 
system with a gapless surface. In the z^ = I case, we can 
of course refermionize the above Hamiltonian and repeat 
the analysis of Sec. IV B to show this explicitly. In the 
fractional case, where refermionization does not result in 
a solvable model, the situation is more subtle as the var¬ 
ious inter-wire terms are irrelevant in the weak coupling 
limit. To avoid additional complications that arise from 
that, we work in the regime where the bare amplitudes 
are large enough such that the inter-wire terms flow to 
the strong coupling limit and successfully gap out the 
bulk. 

Indeed, if the Hamiltonians Hy and Hz have very large 
bare amplitudes, we can neglect the intra-wire terms 
and the physics becomes practically independent of m 
because the inter-wire Hamiltonian is quadratic in the 


fermionic TZ-fields. Since we know from the fermionic lan¬ 
guage that the system is gapped in the m = I case, it is 
clear that in this limit the same is true for m > I. This 
result is expected to remain true for moderately large 
bare amplitudes as well. 

In what follows, the topological nature of the Hamil¬ 
tonian is revealed again with the aid of diagrammatic 
representations. We depict the 77 modes which are not 
gapped by intra-layer terms (i.e., the modes near kx = 0 ) 
by diagrams of the form shown in Fig. (10). As before, 
the symbol 0 ( 0 ) represents a right (left) mover, and 
colored arrows connecting two modes represent coupling 
between them. The horizontal axis represents the layer 
and spin indices (77, s), and the vertical axis represents 
the intra-layer position ( 7 , a). 

Forgetting for a moment that the symbols represent 
dispersing ID modes, if they are treated as states lo¬ 
calized on the corresponding lattice points, the diagram 
presents the 3D problem as a 2D tight-binding model, 
defined in real space by Eqs. (42)-(43). The 2D lattice 
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It It 2t 2t 3t 3t 



FIG. 10. A diagrammatic representation of the reduced 2D 
tight binding model, describing the inter-wire part of the full 
3D model. Here, each “lattice point” corresponds to a lin¬ 
early dispersing mode in the x direction, with a well dehned 
y and z coordinates. The horizontal axis describes the layer 
and spin indices (n, s), and the vertical axis describes the po¬ 
sition within each layer Again, the symbol (Q) 

represents a right (left) mover. The green arrows represent 
the terms in the intra-layer Hamiltonian Hy^ and the brown 
arrows represent inter-layer terms which couple the spin-up 
modes, and take the form —2tz cos (kzo) Tx in /c-space. We 
emphasize that the reduced model dehned in Eq. (42)-(43) 
contains additional terms which are not depicted here. If all 
the terms are considered, the full 3D model is topologically 
non-trivial if the reduced model forms a 2D topological in¬ 
sulator. Thus, the analysis of some aspects of the strongly 
interacting 3D phase is reduced to the analysis of the 2D 
non-interacting topological phase. 


model describing the inter-wire coupling terms is referred 
to as the reduced 2D model. The usefulness of the re¬ 
duced tight binding model description is revealed by not¬ 
ing that the full 3D model is topologically non-trivial if 
the reduced model forms a 2D topological insulator. 

If the system is inhnite (or periodic) in the y and z 
directions, the corresponding momenta are good quan¬ 
tum numbers. We therefore write the problem in terms 
of its Fourier components ky and kz. We note that in 
order to take advantage of the bosonization description, 
we do not perform a Fourier transform in the x direction. 
Consequently, the Hamiltonian takes the form 

H = H, + J2 f (hy (k) + h, (k)) §k(a:), (44) 

k 

with k = {ky, kz), and 

^k(a:) = Y. (45) 


hy (k) = -tyTx - h (jyay + x) (46) 

- y ijx<yx - Tydy) COS i^kyU) 

- ^ {Tx<yy + Tydx) sin {4kya ), 

(k) = 1 [(m' - 2t^ cos (k^a)) (47) 

+ 2t'^ sin (k^a) SyT^] (1 - a^), (48) 

where Ny is the number of wires in each layer, and Nz is 
the number of layers (cf. Eqs. (31) and (32)). 

If the 3D model is finite in the y or z direction, the 
reduced 2D model has an edge. For m' < 2tz, it can be 
checked that the reduced model forms a 2D topological 
insulator, and has counter propagating edge modes. Re¬ 
calling that each lattice point in the reduced model rep¬ 
resents a linearly dispersing r^-mode, we argue that these 
edge states correspond to the gapless surface modes of 
the original 3D model. Therefore, we can get the surface 
modes by diagonalizing the reduced tight binding model 
and finding the corresponding edge states. 

In order to get an explicit analytic form, we focus again 
on the regime where 2tz is close to m', and the gap be¬ 
comes small. If the system terminates at the z = 0 plane, 
the reduced model has two edge modes which take the 
form 

il3{x,ky) = —(49) 
V y j^n 

where /3 = 1,2 labels the two counter-propagating edge 
modes, f{z) is the function defined after Eq. (38), and 
the vectors Ap are 

Ai = i(l + i ,1 — ^ ,—i — l ,i —1 , 

-i + 1 ,-i-l ,-1 + i ,1 + 0 

A2 = AI (50) 


In the integer case, where m = 1, Eq. (49) can be 
thought of as a bosonic description of the Dirac mode lo¬ 
calized near the z = 0 surface (we point out the similarity 
between the fields defined in Eqs. (49) and (38), up to 
a different choice of basis). In the fractional case, where 
m > 1, the above gapless modes cannot be described by 
Dirac’s theory of free fermions. As before, we refer to 
these more general modes as fractional Dirac modes. 

The part of the inter-wire Hamiltonian describing the 
edge modes of the reduced model (which correspond to 
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the surface mode of the original 3D model) takes the form 
Hedge ‘2ityCl ^ ^ dxky^ (x, ky^T^^ix^ ^y) 

ky 

=-2tyi J dxdy^{x,y)fzdy^{x,y), (51) 

where I* = ( ^2 )^, and is a Pauli matrix acting on 

this basis. To simplify the analysis presented in the next 
section, the second line of Eq. (51) was written in the 
continuum limit. 

In the next section we will use the reduced 2D model 
formalism to study the properties of the surface once it 
is gapped by breaking its protecting symmetries. We will 
see that the resulting gapped fractional Dirac modes have 
unique properties which distinguish them from massive 
Dirac fermions. 


D. Gapping the surface 

1. Breaking time-reversal symmetry: Halved fraetional 
quantum Hall effeet on the surfaee 

Having written a low energy effective surface Hamil¬ 
tonian using the reduced model, we turn to study what 
happens when the surface is gapped by breaking time 
reversal symmetry. 

Like in a non-interacting strong topological insulator, 
we expect the effective action describing the response of 
our system to the electromagnetic field to contain a 0- 
term of the form / d^xdtE • B. In general, time 

reversal symmetry severely restricts the possible values 
that the axion angle 0 can take. Naively, one would ex¬ 
pect 0 to be either 0, corresponding to a trivial insula¬ 
tor, or TT, corresponding to a topological insulator [48]. 
However, it was shown in Ref. [18] that fractional axion 
angles may be consistent with time-reversal symmetry in 
topologically ordered systems. As we will see, our frac¬ 
tional strong topological insulator provides an example 
of such a scenario. 

In general, the ^-term in the effective action implies 
that breaking time-reversal symmetry on the surface re¬ 
sults in a non-zero surface Hall conductance of the form 
(Jxy = In what follows, we calculate a^y directly 

and then use the above general result to determine the 
axion angle 0 characterizing our system. To do so, we 
first break time-reversal symmetry by introducing a Zee- 
man field B on the surface. We examine a configuration 
where B changes sign as we cross the line ^ = 0. By 
studying the properties of the gapless mode attached to 
the boundary, we will be able to deduce the surface Hall 
conductance. 

Within the reduced 2D model formalism, the problem 
of finding the ID channel attached to a magnetic domain 
wall on the surface is converted into that of finding the 


localized zero-energy mode associated with a similar do¬ 
main wall on the edge of a 2D topological insulator. To 
be concrete, we use the continuum model described by 
Eq. (51) and add a space dependent perturbation of the 
form where B{y) = Hosign(^). This Hamil¬ 

tonian, and hence the full inter-wire Hamiltonian, has a 
zero-energy solution of the form 

isix) = \/\^a Y^f{n-a) ■ ^j,n{x), (52) 

with 

^=1(010-110-10), (53) 

and As = 2^. 

Notice that this mode, being a combination of right 
movers, is a right moving mode. Eurthermore, we argue 
that it is identical to the chiral mode that resides on the 
edge of a = I/m Laughlin QHE state. 

To see this, we calculate the electron propagator char¬ 
acterizing this ID channel, defined as Gb{x — x'^t — f) = 

. Recall that ^b{x) is an exact zero- 

energy solution of the Inter-wire Hamiltonian for any x. 
Therefore, neglecting the coupling to the gapped excita¬ 
tions of the inter-wire Hamiltonian, we can calculate the 
expectation value with respect to the quadratic intra¬ 
wire Hamiltonian shown in Eq. (41) . This results in the 
propagator 

which is indeed identical the one characterizing the edge 
of a z/ = I/m QHE edge state [49]. 

Invoking the bulk edge correspondence, we find that 
axy = in agreement with Sec. HI. The halved 

fractional surface quantum Hall effect indicates that the 
system is characterized by a fractional axion angle 0 = 
7r/m. A direct implication of the above is that the system 
has a fractional magneto-electric response. 

We note that to get to Eq. (54) we have used 
the full intra-wire Hamiltonian, which includes contri¬ 
butions from modes which are already gapped by the 
inter-wire terms. In doing so, we implicitly assume 
that the intra-wire Hamiltonian does not couple ^b to 
these gapped fields. This can be justified by calculat¬ 
ing ^ (f') , where is one of the gapped modes 

and the average is taken with respect to the intra-wire 
Hamiltonian. Using the definitions of and [Eqs. 
(41) and (52)] and the fact that must be orthogonal to 
a straightforward calculation shows that this average 
indeed vanishes, justifying the above calculation. 
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2. Coupling the surface to a superconductor: The 
emergence of a fractional Majorana mode 


Next, we turn to ask what happens when the surface 
is gapped by breaking charge conservation. This is done 
by proximity coupling the surface to an s-wave supercon¬ 
ductor. It was found in Ref. [39] that when the surface 
of a strong topological insulator is gapped in this fashion, 
the resulting phase resembles a spinless Px + ipy super¬ 
conductor, but has time-reversal symmetry. 

In our model, this problem corresponds to understand¬ 
ing what happens to the edge of the reduced 2D topo¬ 
logical insulator when it is coupled to an s-wave su¬ 
perconductor. This motivates us to write a proximity 

term of the form A ^^1^2 + h.c^^ coupling the two edge 
modes of the reduced model. Using the enlarged basis, 

^ 2 {x,y) ^l{x,y) ) , we rewrite 

the edge Hamiltonian in the form 


1 

2 



-2ityTzazdy H- 


(55) 


In order to reveal the topological nature of the super¬ 
conducting phase, we consider the boundary between a 
region with a non-zero B a region with a non-zero A. 
Studying non-interacting strong topological insulators, 
the authors in Ref. [39] have found that such a boundary 
contains a chiral Majorana mode. We now examine what 
happens in the fractional case. 

In the presence of time-reversal breaking and supercon¬ 
ducting terms, the edge part of the reduced Hamiltonian 
takes the form 


^ J dxdy^;^[-2ityT!,cr!,dy (56) 

C- 


I ^(y)_ ^ I ^(y)_ ^ 

H-H- TzCTx 

a a 


We consider a simple situation, where B{y) = Bo&{y) 
and A = Ao0(—y) (to be concrete, we assume Aq > 
0,Bo > 0). 

The Hamiltonian given by Eq. (56) has a zero energy 
solution of the form 


I 

71 




(57) 


with 


^BA = 


2XbXa 
Xb + Aa 


• a)^ • Tj,n(x)^(ja), 

j,n 

(58) 


9{y) 


e ^^ > 0 
gAA ?/,^<0 ’ 


(59) 


A A = 2 ^^ Xb = 2 ^^- For m = I, we find a self- 
Hermitian combination of right moving fermions, making 
it a chiral Majorana mode, in agreement with Ref. [39]. 
If m > I, we find again a self-Hermitian chiral mode. 
However, in this case it cannot be described by a free 
Majorana theory. In particular, the tunneling density 
of states associated with this mode is proportional to 
as opposed to the free Majorana case, where the 
density of states is constant. This is in fact similar to 
the tunneling density of states characterizing the edge of 
a z/ = I/m Laughlin QHE states. To find this result, one 
can repeat the process that led to Eq. (54) and write the 
propagator of the above fractional Majorana mode. The 
above results are in agreement with Sec. HI, where we 
have used a modified time-reversal symmetry to directly 
model the surface. 


V. DISCUSSION 

The effects of strong interactions on topological matter 
is a complex subject, and a very active field of study. The 
challenge becomes even larger in 3D due to the limited 
analytical tools, and the technical difficulties in apply¬ 
ing the existing numerical methods. The coupled-wires 
approach has been shown to be useful in the theoretical 
study of such phases in 2D. Within this approach, the 
existing machinery for treating ID interacting systems is 
used to produce various analytical results, which reflect 
the topological nature of the phase in question. 

In this paper, we have demonstrated that the coupled- 
wires approach can be used to model and analyze 3D 
topological systems. We have focused on the fractional 
counterpart of the well known strong topological insula¬ 
tor, which demonstrates the remarkable physical prop¬ 
erties of 3D topologically ordered phases. While it is 
not trivial to study the bulk excitations and the gapless 
surface in the current formulation, we were able to an¬ 
alyze the non-trivial characteristics of the surface once 
it is gapped. In particular, the coupled wires approach 
has been shown to be very effective in revealing the na¬ 
ture of ID modes residing in the vicinity of domain walls 
between distinct gapped regions. 

This allowed us to show that if the surface is gapped by 
breaking time-reversal symmetry, it is characterized by a 
halved fractional Hall conductivity of the form (Jxy = 
Eurthermore, if the surface is partitioned into a 
magnetic region and a superconducting region, a novel 
fractional Majorana mode was shown to emerge on the 
boundary between them. 

It would additionally be interesting to generalize the 
configuration discussed in Ref. [50], in which one can 
study and electrically detect the interferometry of Majo¬ 
rana fermions on the surface of a strong topological in¬ 
sulator, to our fractional case. Such an electronic Mach- 
Zehnder interferometer is expected to contain signatures 
that differentiate the fractional Majorana modes from 
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free Majorana modes. This will be elaborated on else¬ 
where. 

The current approach provides a path for construct¬ 
ing the fractional analogues of non-interacting topologi¬ 
cal phases. It would be interesting to systematically ap¬ 
ply this formulation to the different classes found in the 
periodic table of topological phases [13, 14], as was done 
for 2D systems in Ref. [30]. The current work provides 
the necessary tools for extensions to 3D phases. 

Note added. We point out a recent related work (Ref. 
[51]), in which an extension of the coupled-wires ap¬ 
proach to 3D was discussed. We note that the phases 
constructed in the above paper are composed of quan¬ 
tum Hall subsystems, and are therefore fundamentally 
different from the 3D fractional phases constructed in 


the current work. 
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